unclassified 


NAVAL  INIE.LLI6ENCE  SUPPORT  CENTER  WASHINGTON  D C TRA--ETC  F/6  20/4 
LOW-ASPECT-RATIO  WING  IN  A BOUNDED  INVISCID  FLOW  (KRYLO  MALOGO  — ETC(U) 
AU6  76  V I KHOLYAVKOt  Y F USIK 

NISC-TRANS-3838  NL 


I -77 


AD  A 0 32529 


DEPARTMENT  OF  THE  NAVY 
NAVAL  INTELLIGENCE  SUPPORT  CENTER 
TRANSLATION  DIVISION 
4301  SUITLAND  ROAD 
WASHINGTON,  D.C.  20390 


SIFICATION : 


UNCLASSIFIED 


APPROVED  FOR  PUBLIC  RELEASE,  DISTRIBUTION  UNLIMITED 


Low-Aspect-Ratio 


Vling  in  a Bounded  Inviscid  Flow  L - 

Krylo  gaiogo  udlineniya  v .ogranichehnon  j>otoke  nevyazkoy  zhidkosti 


AUTHOR (S) : Kholyavko,  V.I.  and  Usik 


PAGES:  ** 


SOURCE : / Saaoleto3t 
Khar'kov  S 
PP.  3-11 


ORIGINAL  LANGUAGE:  Russia^ 


TRANSLATOR:  C 


N1SC  TRANSLATION  NO 


LOW-ASPECT-RATIO  WING  IN  A BOUNDED  INVISCID  FLOW 


[Kholyavko,  V.  I.  and  Yu.  F.  Usik,  Krylo  malogo  udlineniya  v ogranichennom 
potoke  nevyazkoy  zhidkosti,  in:  Aircraft  Construction  and  Technology  of  the 
Air  Fleet  (Samoletostroyeniye  i Tekhnika  Vozdushnogo  Flota) , No.  20,  Khar'kov 
State  University  Publishing  House,  Khar'kov,  1970,  pp.  3-11;  Russian] 


The  aerodynamic  characteristics  of  a thin  planar  wing  moving  near  a solid  /3* 
or  free  surface  (hydrofoil)  are  determined.  Use  is  made  of  the  thin  body  theory, 
whereby  the  three-dimensional  flow  past  a thin  body  elongated  in  the  direction  of 
motion  is  approximately  replaced  by  two-dimensional  flow  in  transverse  planes. 

In  the  case  under  consideration,  the  problem  is  reduced  to  the  study  of  the  motion 
of  a flat  plate  near  an  interface 


FIGURE  1 


Let  a planar  wing  of  low  aspect  ratio,  whose  maximum  span  coincides  with 
the  trailing  edge,  move  at  a small  angle  of  attack  and  a constant  velocity  V . 

A portion  of  the  wing  dx,  « V^dt  will  pass  through  a fixed  transverse  plane  ZOY 

in  time  dt.  In  the  equivalent  two-dimensional  problem,  this  motion  of  the  wing 
corresponds  to  a vertical  displacement  of  a flat  plate  moving  at  constant  velo- 
city Vq  =-  \'JX  by  an  amount  V^dt  and  to  a change  in  the  width  of  the  plate  from 
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where  £(x^)  is  the  local  wing  semispan. 

The  flow  in  the  ZOY  plane,  caused  by  the  vertical  displacement  of  the  plate 
and  by  the  change  in  its  width,  will  cause  the  apparent  mass  of  the  plate  to 
change  by  an  amount  corresponding  to  the  wing  lift  increment  over  the  length 
dx^,  i.e., 

dY  d (mVo)  _ 1/  dm_(x,) T/  dm  (x,)dx,  _ 1/7  dm  (xj) 

. dxx  dt  dt  - df  - v-a~dir  • 

Integration  of  this  expression  over  the  chord  length  0 < x ^ b gives  the 
wing  lift  e 


Y~y£ldxx=.VlaJn(b), 

0 


(1) 


where  m(b)  is  the  apparent  mass  of  the  plate,  determined  in  the  wing 
section  along  the  maximum  span. 

The  transverse  aerodynamic  moment  relative  to  the  OZ  axis  passing  through 
the  wing  apex  is  calculated  from  the  formula 


- — jx.grfx,  = — Via  |'x1dm(x1)  = — V««x 

o 1 o 

» 

X [m(b)b  — \m(xl)dziy 


(2) 


The  pressure  drag  force  (induced  drag  force)  of  the  wing,  allowing  for 
the  suction  at  the  leading  edges,  is  determined  from  the  relation 


x,  = {y*. 


(3)v 


Instead  of  the  forces  and  moment  (1)  — (3) , we  will  consider  the  following 
coefficients: 


Ct=^  ' .ri 

r X‘  Mx''\ 

c*i  — r-o;  mz.  — 

' 1 <imSbri 


where  S is  the  wing  area  in  the  plane; 


is  the  velocity  liead; 


S = 2| /(*,)<£*,; 


• <7-  " -jp 


Then  the  aerodynamic  characteristics  or  a wing  of  low  aspect  ratio  may  be 
represented  as  follows: 


m. 


2 *[ 1 - tohn  j m • 
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If  the  lift  and  rolling  moment  coefficients  are  known,  the  position  of 
the  wing  pressure  center  relative  to  the  wing  apex  in  fractions  of  the  root 
chord  can  be  calculated  from  the  formula 

WTS}j'nlx,)djc*'  (5) 


xp~  c'  ~ 1 


Formulas  (4)  and  (5)  were  obtained  from  general  relationships  of  the  thin 
body  theory  and  are  applicable  in  cases  where  the  assumptions  of  this  theory  are 
valid.  We  will  note  two  distinctive  features  resulting  from  formulas  (4)  and  (5): 


1.  According  to  Eq.  (4),  the  lift  and  moment  coefficients  are  linear  func- 
tions of  the  angle  of  attack.  The  results  of  experimental  studies  confirm  the 
linear  nature  of  these  relationships  for  wings  of  low  aspect  ratio  up  to  a 4 6°. 


2.  The  lift  coefficient  is  independent  of  the  shape  of  the  wing  in  the 
plane  z^  = £,(x^)  and  is  determined  solely  by  the  wing  cross  section  on  the 

trailing  edge.  Experimental  data  and  calculations  using  exact  theories  also 
confirm  this  characteristic  for  wings  with  aspect  ratio  A 4 1.5. 

Thus,  despite  the  limitation  to  the  case  of  X ->  0,  the  thin  body  theory  /5 

leads  to  qualitatively  accurate  results  for  wings  of  finite  aspect  ratio.  It 
may  be  postulated  that  under  certain  conditions,  this  theory  will  also  yield 
satisfactory  quantitative  results. 

It  follows  from  Eqs.  (3)  and  (4)  that  the  problem  of  determination  of  the 
aerodynamic  characteristics  of  a low-aspect-ratio  wing  in  a bounded  flow  is  re- 
duced to  finding  the  apparent  mass  of  the  plate  near  the  interface.  To  calculate 
the  apparent  mass,  we  will  use  the  method  of  singularities  and  distribute  a 
vortex  layer  of  continuous  strength  y(z)  over  the  length  of  the  plate  within  the 
limits  -l  4 z $ l.  Obviously,  in  this  case  y(z)  = -y(-z) . The  influence  of  the 
interface  will  be  taken  into  account  by  using  the  method  of  specular  reflection 
(Fig.  1). 


The  distribution  y(z)  must  satisfy  the  following  condition  on  the  plate: 
a fluid  particle  adhering  to  the  plate  at  any  point  S must  have  a constant 
vertical  velocity  (Fig.  1).  Using  this  boundary  condition,  we  obtain  an 

integral  equation  for  determining  the  unknown  function  y(z) : 


‘ TU)d2  r 7 <»)(»-*,)  d* 
2-Z,  - (2  - *,)*  + 


2*V0. 


(6) 


Here  the  plus  sign  pertains  to  the  motion  of  the  wing  under  a free  surface 
(hydrofoil),  and  the  minus  sign  pertains  to  the  motion  of  the  wing  near  a solid 
surface  (ground).  The  values  of  £ and  h are  determined  in  the  transverse  flow 
plane  7.0Y  as  a function  of  the  x^  coordinate,  which  enters  into  Eq.  (6)  as  a 

parameter.  To  calculate  the  lift  and  induced  drag  coefficients  of  the  wing,  the 
values  of  £ and  h must  be  determined  in  the  maximum  span  section. 

We  introduce  a new  variable  0 from  the  relation  z = f cosO  (when  -i  4 z £ £., 
we  have  tt  < 0 < 0)  and  make  the  substitution  y(z)  - y(£  cosQ)  - \(0)r  having  set 

li  “ { . Eq.  (G)  is  then  written 
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f 1 (»)  sin  8 d 8 ± p(8)  ( cos  8 — cos  8,)  sin  » _ , 
J cos » — cos 8,  J (Cos  » — cos  «,)»•+  4A>’ 


We  note  the  particular  solution  of  Eqs.  (6)  and  (7)  which  corresponds  to 
the  motion  of  the  wing  in  a boundless  fluid  (h  = °°).  We  obtain  from  Eqs.  (6) 
and  (7) , when  h = ", 


whence* 


f I(£)£2  = f ... = 2icV0 

J J COS  0— cos  d,  O’ 

-/  0 

1 “ 75=5  " 2V"  cl2  8- 


To  solve  Eq.  (7)  in  general  form  for  h / ",  we  represent  the  unknown 
function  y(0)  as  a series 

7(8)  = 2V,,  [j40ctgd+  V A?*  sin  2n8j  . - (9) 


in  which  the  first  term  when  AQ  = 1 is  determined  by  solution  (8),  and  the 

second  term  and  AQ  ^ 1 characterize  the  additional  load  on  the  plate  due  to 

the  influence  of  the  flow  boundary.  It  is  obvious  that  when  h °°,  the  coef- 
ficient Aq  -*■  1,  and  A2n-»-0  (n  = 1,  2...). 

Series  (9)  was  written  by  using  the  condition 
T (»)  = -7  (*-&)- 

If  solution  (9)  is  known,  the  apparent  mass  is  calculated  as  follows. 
The  values  of  the  velocity  potential  function  on  the  surface  of  the  plate  at 
Vq  = 1 are  determined  to  within  a constant  which  later  is  not  significant: 


?u  £-±Tj7(*)*“:Fyj7Wsln»d». 


where  the  upper**  symbol  pertains  to  the  upper  surface  (4>u)  and  the  lower** 
symbol,  to  the  lower  surface  (4’^) • The  general  formula  for  determining  the 
apparent  mass  is^ 


Since  in  this  case,  on  the  lower  surface  of  the  plate  —■  = +1,  and  on  the 
to  dn 

upper  -rf  “ -1,  and  in  addition  dS  = dz,  then 


Translator  s Note:  In  the  display  equations,  read  ctg  as  cotan. 

k 

sic  (translator) 


— — 4p  j <?£  (z)  dz  = — 4p/  j o ^(8)  sin  0 A 0 = 


= 2p  P Jsin8d&$7(&)sin&d0. 


Substituting  the  expression  for  y(0)  from  (9)  into  this  formula,  we  obtain 


m = «p  l*  + y)  • 


Thus,  to  determine  the  apparent  mass  of  the  plate,  it  is  necessary  to  know 
the  first  two  coefficients  AQ  and  A^  of  expansion  (9).  The  remaining  coeffi- 
cients A^n  (n  = 2,  3,  4...)  characterize  the  distribution  of  the  load  over  the 

plate  span  and  do  not  directly  affect  the  overall  aerodynamic  characteristics 
of  the  wing. 

When  h = 00 , it  follows  from  solution  (8)  that  A„  = 1 and  A„  =0,  and 

U 2n 


therefore  by  (10) 


mm  = *p  /*. 


The  aerodynamic  characteristics  of  the  wing,  allowing  for  Eq.  (11),  are 
determined  from  formulas  (4)  and  (5) 

C»--T«  c-i  = yc>-*“icL'  cl 


, * 

XP~  1 — bt * (b)  | ** 


Here  X = — is  the  wing  aspect  ratio, 

= ^(x^)  is  the  equation  of  the  wing  planform. 

Formulas  (12)  are  known  relationships  for  a wing  of  low  aspect  ratio  in  an 
unbounded  fluid. 

Let  us  now  calculate  the  expansion  coefficients  A^  (n  = 0,  1,  2...).  In  /7 
view  of  (9),  Eq.  (7)  is  transformed  as  follows: 

m 

— D <42ncos2rt8f  ± /(it,  ft.)  = 1,  (13) 

fl**l 

wliere  the  following  symbols  have  been  introduced 
/ (ft.  »,)  - A„  (i  Jt  -cos  a, 

yCOsS,  V Ain  + j?n\l)  4'  ~r  ^^2n(^2/i-2  — «^!/i+2); 

n-l  «i=.| 

y _lTcosm8d8  .* 

e J (cos  8 — cos  8,?  -t-  4 n- 

After  calculation  of  integrals  J and  further*  transformations  of  formula 

m 

(13),  we  finally  obtain  equations  for  the  expansion  coefficients  A„  (n  = 0,  1.  2....) 

2n  ’ 
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1.  For  the  motion  of  a wing  near  a solid  surface 


I 


a* 

A0{.0(ii,  8,)  + 1 = 2 A2„[(.2n(Ji,  8,)  — cos 2/i 8,].  (14) 

#X*1 

2.  For  the  motion  of  a wing  above  a free  surface  (hydrofoil) 

At  [2  + C 0(h,  8,)]  — 1 = S Azn  [C2„(ft,  »,)  -cos2/i8,].  (15) 

0=1 

The  following  symbols  have  been  introduced  into  Eqs.  (14)  and  (15): 
C 0(h’  8t)  = — £-|cos8,|  — 2hy\ 

8,)=( — S-J  Tin  (t); 

u = ^Vv  — (sin2  8,  + Ah2} ; 

u = ±VV-+  sin*  8,  -f-  4ft* ; 

V = ]/  (sin!8,  + 4 A*)*  + 16  h*  cos*  8,; 

t = -y-JH— ; Tin  (0  = COS  (2 n arc  cos  /), 

Vu?+~' 


where  T2n(t)  is  a Chebyshev  polynomial  of  the  first  kind. 

If  the  series  of  0g  values  are  fixed  in  relation  (14)  or  (15) , systems 
of  algebraic  equations  are  obtained  for  determining  the  expansion  coefficients 
A2n  = I***)* 

The  number  of  fixed  points  0 ^ (i  = 1,  2...)  determines  the  number  of 

equations  in  systems  (14)  and  (15),  and  hence,  the  number  of  desired  coefficients 
A^n  (n  = 0,  1,  2,...,  i - 1).  The  remaining  coefficients  (n  = i;  i + 1,...) 

must  be  taken  to  be  equal  to  zero. 

* Equations  (14)  and  (15)  were  used  to  calculate  the  coefficients  when  a 
single  point  (approximation  I),  three  points  (approximation  II)  and  five  points 
(approximation  III)  are  fixed. 

Part  of  the  calculations  for  the  motion  of  a wing  near  a solid  surface  /8 
are  shown  in  Fig.  2.  From  known  coefficients  and  A^  and  formula  (10),  one 

can  find  the  apparent  mass  of  the  plate,  and  from  formulas  (4)  and  (5),  the 
aerodynamic  characteristics  of  the  wing. 

Figure  3 shows  values  of  the  derivative  of  the  lift  coefficient  with  respect 
to  the  angle  of  attack  Ca , referred  to  the  value  of  Cn  in  an  unbounded  flow 

y _a  y« 

(12).  Points  denote  average  values  C , obtained  from  the  theory  of  a vortical 

y 4 

lifting  surface  for  rectangular  wings  with  X - 0.4-1. 6. 

Analysis  of  the  calculations  performed  shows  that  satisfactory  results  up 

to  h 0.1  can  already  be  obtained  in  approximation  III  for  a wing  near  a solid 
surface  and  in  II  for  a hydrofoil. 
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FIGURE  3 


If  the  treatment  is  limited  to  the  values  h £-0.4  (wing  near  the  ground)  and 

h £ 0.25  (hydrofoil),  the  aerodynamic  characteristics  of  a wing  of  low-aspect- 
ratio  can  be  determined  from  approximation  I.  The  simple  analytical  relations 

which  follow  from  (14)  and  (15)  hold  in  this  case  when  0 = ? and  A„  =0 

(n  = 1,  2,  3...).  s 2 2n 

1.  For  a wing  near  a solid  surface 


_V  \ + aK':  _ „Vl+4h* 

A6  = — , m » rp/J  - - ? . 


2/i 


Since 


^ _ Vl  + 4A*' 


2.  For  a hydrofoil 


• ; m — vpP  ■ 


V\  + 4h* 


VTT^  t 

Calculations  made  by  using  formulas  (16)  and  (17)  are  indicated  by  solid  /9 
lines  in  Fig.  2. 

Let  us  note  that  if  2 is  replaced  by  a higher  number  (for  example,  3)  in 


the  expression  — 1 , formulas  (16)  and  (17) 

2h 


give  satisfactory  agree- 


ment with  the  calculations  made  in  approximation  III  up  to  h >/  0.10. 

Thus,  the  following  approximate  formulas  can  be  recommended  for  calculating 
the  apparent  mass  of  a plate  in  a bounded  flow  for  h > 0.10: 


1.  Near  a solid  surface 


Vl  + 9h*  n „ox 

m = 3S (8 


2.  Under  a free  surface 


V i -e  9A* 

Formulas  (18)  and  (19)  are  convenient  for  calculating  the  moment  character- 
istics and  position  of  wing  pressure  centers.  For  example,  for  a wing  near  a 
solid  surface,  we  find  from  (5)  and  (18) 


< P=  1 P V r~<Xl)  + 94—  dxr- 

V bl*(b)V  1 + 9A*  J » 


Outside  the  integral  h - h is  the  position  of  the  wing  trailing  edge 

relative  to  the  interfaces,  and  £(b)  is  the  wing  semispan. 


attack 


Jn  the  integrand,  h and  £ depend  on  the  coordinate.  At  small  angles  of 


h (x,)  ~h  + (b-Xi)*~l(b)  [/I  H-  ~ (1  - xja 


8 


Under  the  assumptions  of  the  thin  body  theory 


, and  under  certain 


_b 2. 

» 1(b)  \ 


conditions,  the  second  term  on  the  right-hand  side  may  be  of  the  order  of  the 
first  term;  therefore,  the  position  of  the  pressure  center  on  a planar  wing  near 
the  interface  depends  on  the  angles  of  attack. 


We  will  confine  the  discussion  to  the  case  of  a -*■  0;  then  h(x^)  = £(b)h  = 
const,  and  the  position  of  the  pressure  center  is  determined  from  the  formula 

xn=  1 ” I {Xl)  V Jt  + 9/r*  (20) 


— ■ I (Xl)  — ^1 

where  / (xt)  y is  the  relative  local  wing  semispan;  x^  = — . 

We  will  consider  a family  of  wings  whose  planforms  are  given  by  the  / 

equation  i-(x^)  = x™  . 

The  exponent  m changes  within  the  limits  0 < m < °°.  For  in  = 1,  we  obtain  a 
triangular  wing,  and  for  m -*•  0,  a rectangular  one. 


By  (20),  we  have 


XP=  1 


Vi  + 


1 + 9A2  J 


If  such  transformations  are  made  for  a hydrofoil,  it  follows  from  (5)  and 
(19)  when  a 0 that 

/ 3A  \ nl*(x1)VPCxl)  + 9P  '1 

*p=  1 — \2~y  1+9A2  IJ  2V  /3(ix)+9A*-3  h (2?) 


m 


For  the  family  of  wings  £(x^)  = x^ 

‘ i]m  Vxfm  + 9/7 


*„=  1 — 1 2 — 


3/i 


V 1 +9As  I]  2Vx\m  + 9hi -Zh 


T~dx1. 


(23), 

■■■4 


for  h = 00  (flow  of  unbounded  fluid) , by  (21)  and  (23) 


2m 

Jtp~  ~ 2/n  + 1 ' 


(24) 


This  result  can  also  be  determined  from  formula  (12) . 


In  the  general  case,  the  integrals  in  formulas  (21)  and  (23)  are  expressed 

3 XP 

by  hypergeometric  functions.  Figure  4 shows  the  calculations  of  x = — for 

P Xpoo 

1 * 
wings  with  m = — , 1,  2 (solid  lines  - for  a wing  near  a solid  surface;  dashed 

lines  - for  a hydrofoil) . 


Analysis  of  the  results  obtained  shows  that  as  the  wing  approaches  a solid 
surface  (ground),  the  pressure  center  shifts  toward  the  trailing  edge,  and  this 
shift  is  greater  the  smaller  the  exponent  m.  In  particular,  the  greatest  shift 
should  be  observed  for  a rectangular  wing  (m  ->  0) . 

The  influence  of  the  free  surface  boundary  for  a hydrofoil  is  the  opposite 
of  the  ground  effect. 
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